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Abstract 

Similar to the associated Legendre functions, the differential equation for the associ- 
ated Bessel functions Bi m (x) is introduced so that its form remains invariant under 
the transformation I — > — I — 1. A Rodrigues formula for the associated Bessel func- 
tions as squared integrable solutions in both regions I < and I > is presented. 
The functions with the same m but with different positive and negative values of I are 
not independent of each other, while the functions with the same I + m (I — m) but 
with different values of I and m are independent of each other. So, all the functions 
Bi,m[x) may be taken into account as the union of the increasing (decreasing) infinite 
sequences with respect to I. It is shown that two new different types of exponential 
generating functions are attributed to the associated Bessel functions corresponding to 
these rearranged sequences. 
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1 Introduction and Motivation 

The generating functions have found many applications in the physical, chemical and math- 
ematical systems. The study of the important chance process called the branching process 
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[D 121 E], random graphs and complex networks [1], polymerization kinetics [5], counting 
problems in combinatorics [6], are some applications of the theory of generating functions. 
The generating functions are used to obtain expected values (averages), variances, moments 
and cumulants of distributions, and also, to establish relationships between distributions [7]. 

The exponential generating functions and their numerous generalizations have been al- 
ternatively introduced and studied by various methods for the orthogonal polynomials and 
special functions (see Refs. [H El [lOl HH H2l H31 HH [151 HS1 HZl HH HH [2UI ED [22]). The 
generating functions as continuous functions generally describe the convergence of an infinite 
summation of given infinite sequences of functions. In this sense, the special polynomials 
and functions of a given sequence in one variable x can be defined as the coefficients in 
the expansion of their generating functions. The generating functions include various useful 
properties and all information that is needed to generate the solutions corresponding to a 
differential equation or a set of recursion relations between those solutions. Therefore, gen- 
erating functions are very useful to analyze problems involving summations on the infinite 
sequences of functions such as coherent states. The application of generating functions for 
known orthonormal special functions, allows one to derive a compact formula for the co- 
herent states. Generating function corresponding to a given set of special functions is not 
unique. This manuscript has been devoted to introducing new generating functions for the 
solutions of the differential equation of associated Bessel functions which can be applied to 
obtain bound states of some solvable models in the framework of supersymmetric quantum 
mechanics, such as radial bound states of the hydrogen-like atoms [23] ■ Then, let us re- 
member that Krall and Frink have first studied the Bessel polynomials in the formalism of 
hypergeometric functions [21]. Also, some authors have introduced some generating func- 
tions for Bessel polynomials [23 [26j [27J • Moreover, the generating functions associated with 
the group-theoretic techniques and the Stirling numbers of the second kind have been derived 
for the special sequences of the generalized Bessel polynomials [2S1 [2S1 ED] • Two different 
types of g-analogues of the generating functions for generalized Bessel polynomials have been 
calculated in Ref. [31] too. 

In comparison with the Bessel and Romanovski differential equations [32], the Hermite, 
Laguerre and Jacobi ones from the viewpoint of their polynomial solutions application to 
the physics problems, have attracted much attention until now. However, the Bessel and 
Romanovski polynomials have also been applied to obtain the wavefunctions of some of the 
physical potentials. For example we mention the factorization methods for the differential 
equation of associated Bessel polynomials which are applied to obtain the supersymmet- 
ric structures corresponding to the radial bound states of the hydrogen atom [23]. The 
solutions of the Schrodinger equation for some noncenteral potentials such as hyperbolic 
Scarf and trigonometric Rosen- Morse [3"2"l [331 [31], and of the Klein-Gordon equation with 
scalar and vector potentials, are obtained in terms of the Romanovski polynomials |35j . The 
trigonometric Rosen- Morse potential is an appropriate candidate to describe the quark-gluon 
dynamics in QCD, since, it can be considered as an appropriate approximation of the combi- 
nation of the Coulomb, the infinite wall and the linear potentials. Therefore, the Bessel and 
Romanovski nonclassical polynomials have the merit of taking more into account in order to 
derive new relations. 

The manuscript has 2 sections. In Section 2, we introduce the differential equation of 
associated Bessel functions B[ q ^\x) in terms of the indices / and m, in similarity with 
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associated Legendre functions Pi >m (x): 



(l-x 2 )P;[ m (x)-2xPi m (x) + 

Despite of the Pj m (x)'s, the differential equation for the associated Bessel functions is altered 
when m is replaced by — m. In conclusion, the function B^f^x) is not another solution for 
it. However, the differential equation is invariant under transformation I — > — I — 1, and it, 
in turn, leads to considering the solutions as B\ q ^\x) with / < 0, too. Consequently, we can 

offer a Rodrigues formula for B[ q ^\x) as the squared integrable solutions with I < and 
/ > 0. Furthermore, simultaneous realization of laddering equations with respect to I and m 
by the given Rodrigues formula in both regions, is considered. In Section 3, it is shown that 
the independent solutions B\ q ^ (x) are classified in three different types of infinite sequences. 
The first type sequences of the associated Bessel functions have the same / but different m. 
The generating functions corresponding to them can be followed via the known generating 
functions of the Bessel polynomials. The second and the third types of the sequences are 
constituted by the independent associated Bessel functions with the same / + m and / — m, 
respectively. These two later types provide the possibility to view the set of the independent 
associated Bessel functions in two new perspectives, different from the first one. Finally, we 
calculate two new kinds of generating functions for any of these types of sequences depending 
on whether / + m and I — m are odd or even. 



1(1 + 1)- 



m 



Pi, 



0. 



2 Associated Bessel functions 

This section includes, in addition to review some results of Ref. [36], the formulation of 
Rodrigues representation for the associated Bessel functions B\ q ^\x) with I < 0. Let us 
first remind that the generalized Bessel polynomials of degree n [37J, i.e. 

B^(x) = (jX (> +2 ^) = a n (a,^y n (x;a + 2,P), (1) 



x"e * 

are eigenfunctions of the following linear second order differential operator: 

X ~ a ^Tx ( X " +2e ^) B n^\x)=<n + a + l)B^\x), (2) 

where a n (a, /3)'s are the normalization coefficients. It must be pointed out that the y n (x; a, (3) 
representation of the generalized Bessel polynomials is given by [29, [30] 



»- 1 ( l ) ( ° + " t" ~ 2 ) «(!) k -(f ) " (!)■<» 

For (3 > and a < —2, the generalized Bessel polynomials are orthoghonal and square 
integrable with respect to the weight function x a e^ in the interval < x < oo. Choosing 
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n = I — m + | and a = 2m — q in the differential equation (2), it is straightforward to show 
that the associated Bessel functions [361 



Bl > m {X) a,_ m+f (2m-^) X ^- + f (X) - UJ ^ ) ' ( } 

with / — m + |>0asa non-negative integer, satisfy the following differential equation 
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(9,0) 



C j = 0. (5) 



Similar to the associated Legendre differential equation, since equation (5) is unaltered when 
I is replaced by — I — 1, the functions B\ q ^ (x) with negative / are also another solutions for it. 
But contrary to the associated Legendre differential equation, which is quadratic in terms of 
m, the associated Bessel differential equation (5) is linear in terms of it. It will be clear from 
our discussions that the functions B^' with m < I + | and m > — I + 1 are not normalized 

by the weight function x~ q e^~. Therefore, we limit our study torn— |</<| — m — 1 with 
m < in which q is an integer number. In Fig. 1, we have schematically shown the 2- fold 
hierarchy of all the associated Bessel functions corresponding to q = 6 as the points (l,m) 
on a flat plane whose horizontal and vertical axes are labeled with / and m, respectively. 
Note that for an odd integer q, one of the two parameters / and m has to be half-integer, 
and the other has to be integer. The Rodrigues formula (4) shows that the associated Bessel 
functions are finite summations of (not necessarily positive) integer and half-integer 

powers of x when q is even and odd, respectively. 

Now, in order to formulate the laddering relations with respect to I and m, and also to 
realize the square integrability condition, it is necessary that we obtain the highest powers 
of x in the associated Bessel functions B\ q ^\x) for / < and / > 0, respectively: 

' a t , m ( g ,^)(-l)- m - t+ ^V a+i r( rgtlr^ "^^ 1 + 0(x- 1 ^- 2 ) l>0. (6b) 



Note that (6a) is directly derived using (4), while (6b) is followed from the following relation 
and also (6a), 



d_ 

dx 



21+1 

21 -£ 

x e * 



) = /3 2 '+ 1 x- 2 '- 2 e"^ / > 0. (7) 



The inductive reasoning can be applied to prove the relation (7). The orthoghonality of the 
associated Bessel functions for a given m, and also their square integrability for / < and 
/ > are obtained, respectively, as 

Btf'B^A, = *, a Uq , ^) rH 7_:j ) g- T m i) +f + 1) . (8) 

They follow from the relations (4), (6) and applying integration by parts / — m + | and 
— I — m + | — 1 times, respectively. In the /c-th stage of these processes, the total differential 
expressions become zero, because V + m — | + + 1 < and —V + m — | + A;<0. 
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Similar to Ref. [36J, the associated Bessel differential equation (5) can be simultaneously 
factorized by the ladder operators 



i± -±x* d i (iT q )x± i±imM 

%™~ ±x dx + Y T 2J X± 21 



.4- d (3 8 



? 7 \ /i ,9 



with the following eigenvalues 

(l-m+%)(-l-m + $) B 2 
Ehm = - 2j \ 2 £ l , m =[±-l-m){l-rn+± + l), (10) 

as shape invariance symmetry equations for the indices (l,m) and (I — l,m) as well as (l,m) 
and (l,m — 1), respectively: 

Kn*; m B&%) = E hm B%p(z) A^A+ m B^ m (x) = E, m B^Jx), (11a) 

KA^\x) = £ l>m B^(x) A'A+B^ix) = S l>m Bi%^(x). (lib) 

In this paper, we are interested to present the raising and lowering relations of the indices I 
and m of the associated Bessel functions for both regions, I < and / >0ofm — |</< 



| — m — 1. The shape invariance equations with respect to / and m are realized for every 
normalization coefficient. However, realization of the laddering equations 

At, m BtZ(x) = y/E^B%P{x) A- m B^\x) = y/E^B^x), (12a) 

AiBl%%(x) = y^B^\x) A-B^\x) = y/S^B^x), (12b) 

imposes two recursion relations with respect to I and m, respectively, on the coefficients 
a i,m(q-> 0). One can easily show that the laddering equations (12a) and (12b) as well as their 
corresponding recursion relations for the normalization coefficients ai }Tn (q,/3), are simultane- 
ously established if the latter are chosen as 



^ i( -\ )S ; m = m-|</<0 (13a) 

-3+i)r(-i-m+$) 2 - 

/ i \ — I— m+% — 1 

[ — : 0<Z<|-m-l, (13b) 



lv( ^-m+2+l N )^('-^-m^2^ 



2 J \ 2 , 



where m < One cannot actually follow the raising and lowering form from shape 

invariance symmetries by means of the arbitrary normalization coefficients. The point is 
in order to realize the laddering symmetries (12a) and (12b), it is a necessary condition to 
select the normalization coefficients as in (13a) and (13b). 

Our reason for making the associated Bessel functions with I < goes back to the fact that 
the functions Bq^\x) with m < ^y- are not annihilated by the operators AQ m . Therefore, 
for a given m, decreasing of the index I can be terminated at I = m— |. Indeed, from the Eqs. 

(12a) we have At Bq 1 '^ , (x) = and A~ q B^' 13 } (x) = 0. Moreover, from the Eqs. 

(12b), it becomes clear that the associated Bessel functions lain on the lines I = m — | and 
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/ = — m + f — 1 are annihilated by the ladder operators shifting m: A^ + q +l B^ q '_^l(x) = and 
^i-i-^ii-t-i ( x ) = 0- According to the above discussions, for a given nonnegative I and for 

2 '2 

every m < we have B_l_ lm (x) = /?( — B\ q ^\x) , which means the functions with 
/ < and / > settled on the horizontal lines of the Fig. 1 are mutually dependent on each 
other. However, each of the vertical and oblique lines possess the associated Bessel functions 
that are independent of each other. This allows us to construct the generating functions 
for them in three different methods by using Rodrigues formula (4) for both regions I < 
and / > 0. Note that the differential equation (5) has two independent solutions and above 
discussions have focused on one of them. 



3 Exponential generating functions for the associated 
Bessel functions 

The square integrable associate Bessel functions can be applied to obtain bound states corre- 
sponding to some one-dimensional supersymmetric potentials and also some two-dimensional 
quantum mechanical models having a Lie algebra symmetry [231 EH] - Therefore, exponen- 
tial generating functions corresponding to the formal power series of associate Bessel func- 
tions Bi 9 J 3 \x) with the same /, the same / + m and the same / — m, are important not 
only from the point of view of mathematical derivation but also from the point of view 
of physical applications. What we do is to consider three different methods for comput- 
ing the generating functions, based on the presentation of appropriate infinite sequences of 
the associated Bessel functions. The first type of the infinite sequences is {Bl q J 3 \x)}^ 1+l 

for / < or {Bi 9 J 3 \x)}^ 1 _ l _ 1 for / > 0. Using the definitions n : = | — I — m — 1 

and p := I — m + ~, all the associated Bessel functions B\ q ^\x) with / < and / > 
can also be rearranged as the union of the second and third type infinite sequences: 
{B < f^_ n+ q_ l {x)} c ^ =Q and {B^f^ +1 {x)}^ =Q ^ respectively. The second type sequences are in- 
creasing with respect to the index / of the associated Bessel functions Bl q n f \x) while the 
third type sequences are decreasing. These sequences are automatically covered by the as- 
sociated Bessel functions which are linearly independent on the interval < x < oo with 
respect to the inner product (8), as {B^^^ Jx)}^^ \J{ B f-i-2k-i,J x )}~rZ^-k-i 
and {B^'fl g ,. (x)}-°° q , . \]{B {q 'f}. q (x)}-°° q . , with k = 0, 1, 2, • • • . Now we are 

L m+2k— j+l,m\ ' J m=^ — fe— 1 L m+2k—^,m^ 'Jm=j— k— 1 

in a position that for the sequences given in above we can derive three different types of the 
generating functions. 

Generating functions for given q and I: Let us first introduce the generating func- 
tions for the power series with different m but with the same /, which are calculated in a way 
similar to those for generalized Bessel polynomials [281 1291 130] . The generating functions 
corresponding to the first type sequences of associated Bessel functions with the different m 
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but the same /, 



for I < (14a) 
for I > 0, (14b) 



for \t\ < 1, are 

d(x, t) = (l + t) 21 x z+ 5 e ^TTi) + li m 

Z— Yx 

= (1 + tf l x l+ ie^) . 

The relation (15), for both regions I < and I > 0, are followed by the following Cauchy's 
integral formula 

tv, 1 =— l dz ,*" 6 * (16) 

r(k + l)\dx) \ J 2mJ c{X)t) (z-x)^ 1 ; 

with k — I — m + |. The Fig. 2 shows that the integration contour C(x, t) is a closed path 
around the circle \z — x\ = \t\x in positive direction, z = has also been settled out of the 
contour C(x, t). Also, the real pole z = x(l + 1) is always settled inside the contour, since it 
is on the circle at the right and left hand sides of center whether t is positive or negative. 
Considering the following relation 

*^(*) = «Ug,fl(-i^ (£) > ( 17 ) 

one can establish the connection between the Gi(x,t) with the Laguerre generating function 
(1.18) of Ref. [30]. 

Generating functions for given q and I + m: In this case the sequences are increas- 
ing with respect to the I. Due to the fact that whether n is odd or even, i.e. n = 2k + 1 or 
n = 2k, highest functions are -B_l_l s._ k _ 1 {x) and B^fjfl^Jx), respectively. These func- 
tions lie on the lines m — I + | and m = I + | — 1 of Fig. 1, respectively. Therefore, 
it is obvious that the terminology of highest functions has been devoted to the associated 
Bessel functions B ( f^ i \x) with the most value for m. a) First we suppose that n is odd, i.e. 
n = 2k + 1. For a given value of k, the generating functions corresponding to the second 



Gi{x,t) :-- 



E 
E 



m=0 m! a 

oo i 
m=l m! a 



J (,,/J) Ax) 



i,/ — m+ 



(<?,/?) 



? (9 ' ,3) ,(x) 



T(2/ + l) 



dz 



21 



z 2l e * 
z- x(t + 1) 



(15) 
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type sequences are calculated as 



G n =2k+l(x, t) — 



B 



(q, 



m— fc— 1,2— m- 



= X 2 



' Q ( 2m V- Om-fc-l,a-m-fc-l(9>/3) 



e 1 



m=0 

oo 



£ 2 



(xt) m 
^ 2ttz 

m=0 



^2m-2fe-2 



e * 



<7(x,t) 



47ri xyxt 
fc-i 



2; 



-2fc-2/ 



(2 — x) 
— x)e" 



2m+l 



1 , 

2 



C(x,t) " l-Vxt 

.\ 2fe+l 

+ 



dz 

C(x,t) 
\ 2fc+l 



z 2fe 2 (z — x)e~ 
z — 



l+Vxi 



(18) 



As it has been shown in Fig. 3, C(x,t) is a closed contour in positive direction around a 
circle with center at (jz^, 0) and radius of R — on the complex plane z. Again, z = 
is out of the contour C(x,t). Also, the variable t is considered as t < -. Therefore, z = x 
settles on the real axis at the inside of the circle and at the left hand side of the center. 
Moreover, the poles z = - V and z = - X r- have been settled on the real axis as well as 

1 1 l — Vxt l+Vxt 

on the circle, at the right and left hand sides of center, respectively, b) Second we suppose 
that n is even, i.e. n = 2k. For a given value of k, once again using the integration contour 
C(x,t), the generating functions of the second type sequences are calculated as 



Gn=2k(x,t) — 



B {q,P L „ Ax) 

m— — m— k— 1 \ / 



" Q (2m + 1)! o m _ fc) s_ ro _ fc _i(g, (3) 

(l - Vxiy k - (l + Vxtj 



S.—k 
X2 



2Vxt 



2k 



(19) 



Note that the series in (18) and (19) are summed on the parameter m for given values 
I — 2k — 2 and | — 2k — 1 of I + m, respectively. 



Generating functions for given q and I — m: This case involves the decreasing se- 
quences with respect to the / and again, there exists two possibilities: p can be both odd 
and even. For p = 2k + 1 and p = 2k, the highest functions corresponding to them are 
B^'H k _A x ) an d B^_f \_ k _A x )> respectively. Also, these functions have been placed on the 
lines m — — I + | — 1 and m = — I + | — 2 of Fig. 1, respectively, a) In the case that p is 
odd, i.e. p — 2k + 1, for a given value of k and for \t\ < 00, the generating functions of third 
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type sequences are calculated as follows 



G p= 2k+l(x,t) — 



oo g^P) | 

t k—m,^—m—k—l 



X 



( ,_ ml a k _ m z_ m _ k _ x {q, ft) 



(2k + l)\x 



2-1-1.4-1 £ 



E 



(xt) 7 



d 

dz 



m=0 
2fc+l 



<"/2 



2m ml J c{x) (z - x) 



Jlk-2m 

z e z 

_ T \2k+2 



z 2k e^ * 



(20) 



In order to satisfy Eq. (20), it is sufficient that the arbitrary contour C(x) is chosen so that 
the points z = x and z = lay inside and outside of that, respectively (see Fig. 4). b) 
When p is an odd number, i.e. p = 2k, the generating functions of third type sequences for 
a given k are calculated as follows 



G p=2 k(x,t) = 



I fe— m— m— fc— lv ' 

( m! a A ._ m _ li |_ m „ fc _ 1 (g, /3) 



(2fc)!x 



E 



(xt) 7 



m=0 ^0(a;; 



^ c 2 



,2fc+l 



X 2 



+fc+l 



e * 



f/2 



2/v 



„2fc-2 



(21) 



Here, the series in (20) and (21) are summed on the parameter m for given values 2k — | + 1 
and 2A; — | of Z — m, respectively. The accordance of the above generating functions with 



the theorem 1 of Ref. [29J can be considered as new confirmation for it. 

If we choose q = 0, then we can claim that the relations (18), (19), (20) and (21) are 
generating functions corresponding to the associated Bessel functions with l+m = — 2(k+l), 
I + m = —2k — 1, I — m = 2k + 1 and I — m = 2k, respectively. Therefore, we have obtained 
four new different types of generating functions for the associated Bessel functions depending 
on whether I + m and I — m are negative even or negative odd integers and positive odd 
or nonnegative even integers, respectively. Therefore, in order to obtain new generating 
functions we have used square integrable associated Bessel functions in both regions I < 
and I > with the same Rodrigues representations for them. 
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Fig. 1. The comprehensive plan of the squared integrable solutions for the differential equa- 
tion (3) of associated Bessel functions with q — 6. 



Fig. 2. Plot of the integration contour C(x, t) for the generating functions of first type 
sequences of the associated Bessel functions. 



Fig. 3. Plot of the integration contour C(x,t) for the generating functions of second type 
sequences of the associated Bessel functions. 



Fig. 4. Plot of the integration contour C(x) for the generating functions of third type 
sequences of the associated Bessel functions. 
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